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‘Report: Carbon CycCiny in 
‘Hip h -Latitude ‘Ecosystems 

Alan Townsend, Stephen Frolking, and Elizabeth Holland 

The carbon-rich soils and peatlands of high-latitude ecosystems 
could substantially influence atmospheric concentrations of CO a 
and CH 4 in a changing climate. Currently, cold, often waterlogged 
conditions retard decomposition, and release of carbon back to the 
atmosphere may be further slowed by physical protection of organic 
matter in permafrost. As a result, many northern ecosystems accu- 
mulate carbon over time (Billings et ah, 1982; Poole and Miller, 
1982), and although such rates of accumulation are low, thousands 
of years of development have left Arctic ecosystems with an 
extremely high soil carbon content; Schlesinger’s (1984) average 
value of 20.4 kg C/m 2 leads to a global estimate of 163 x 10 15 g C. 

All GCM simulations of a doubled C0 2 climate predict the greatest 
warming to occur in the polar regions (Dickinson, 1986; Mitchell 
1989). Given the extensive northern carbon pools and the strong 
sensitivity of decomposition processes to temperature, even a slight 
warming of the soil could dramatically alter the carbon balance of 
Arctic ecosystems. If warming accelerates rates of decomposition 
more than rates of primaiy production, a sizeable additional accu- 
mulation of C0 2 in the atmosphere could occur. Furthermore, CH 4 
produced in anaerobic soils and peatlands of the Arctic already com- 
poses a good percentage of the global efflux (Cicerone and Oremlund, 
1988); if northern soils become warmer and wetter as a whole, CH 4 
emissions could dramatically rise. A robust understanding of the pri- 
mary controls of carbon fluxes in Arctic ecosystems is critical. 

As a framework for a systematic examination of these controls, we 
discussed a conceptual model of regional-scale Arctic carbon 


316 


Modeling the 'EarthSystem 


turnover, including CH 4 production, proposed by E. Holland, and 
based upon an extension of the Century soil organic matter model 
(Parton et al„ 1987, this volume). The details of Century will not be 
repeated here; rather, we will restrict our discussion to the specific 
modeling challenges posed by Arctic ecosystems. 

Biophysical Model 

Both soil temperature and soil moisture status are critical deter- 
minants of carbon dynamics in any ecosystem. Some unique fea- 
tures of the Arctic physical environment complicate the prediction ot 

these variables: 

. Permafrost, or permanently frozen subsurface soil, underlies 
much of the Arctic. It has a significant effect on regional hydrol- 
ogy by providing an impermeable barrier to vertical infiltration. 
Changes in permafrost will alter the size of the active carbon pool, 
regional hydrology, and perhaps topography through uneven set- 
tling The extent to which the permafrost layer may change In a 
warmer climate is not well known, but recent evidence from oil 
wells in northern Alaska suggests a distinct warming trend in the 
permafrost over the 20th century. 

• Snow, being a very effective insulator, can dramatically affect the 
soil thermal regime. In addition, snowmelt is a dominant source 
of water to Arctic ecosystems. 

• Freeze/thaw processes in the soil are Important to the regional 
hydrology (frozen soil is impermeable) and to the soil therma 
regime due to delays at 0°C resulting from the energy require- 
ments of a change in phase. The timing of the spring ^w Is cru- 
cial, because it occurs when available sunlight is high ( ay 
June). An earlier thaw would greatly enhance ecosystem 
productivity, and a deeper thaw in permafrost regions would 
thicken the active decomposition layer. 

Due to high spatial variability of these factors in the Arctic, we 
believe a fairly comprehensive biophysical model of the soil environ- 
ment that provides input to the carbon turnover model will be nec- 
essary for regional-scale simulations (Figure 1). The present concep- 
tual structure of this model contains four distinct layers (Figure 2): 
unsaturated peat or soil at the surface; saturated peat or soil, deep, 
unfrozen, saturated mineral soil or extensively decomposed peat; 
and permafrost. Changes in the extent of the top two layers (that , 
changes in the water table) may occur on time scales of days to 
weeks, while significant changes in permafrost are only likely to 
occur on decadal time scales. 
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fhl o .I*", rr , ‘ u ^ctic carbon cycling model with 

he sod biophysical model. Inputs are climate variables, intermediate results 

e sod temperature and moisture profiles , and final outputs are C0 2 and 
CH 4 fiux time series and the annual change in sod carbon storage. 

Model inputs include the thermal and hydraulic properties of the 
surface vegetation and soil layers, and a parameterization of lateral 
surface water movement. The model will be driven by (at minimum) 
air temperature and precipitation, and will determine the tempera- 
ture profile of each layer, the thickness of the aerobic vs. saturated 
ayers, seasona! depth of thaw, and changes in permafrost. As well 
the model will need to determine snow depth, extent of snow cover 
and timing of snowmelt. 

Carbon Turnover Model 

Of So V H 1 l ° pmei ? t ° f a good biophysical model is just one of a number 
of modeling challenges posed by Arctic ecosystems. While fluxes of 
carbon between the atmosphere and ecosystems are low, the stocks 
o available substrate in the soils are enormous. As a result the ini- 
tial response of Arctic systems to climate change will not be dra- 
matic, but the cumulative effects and feedbacks to the climate sys- 
tem on longer time scales are potentially extreme. The most critical 
unknowns in predicting these dynamics concern the responses of 





heat and 
moisture 
transfers 


Saturated 

peat/soil 


Temp(t) 

thickness(t) 



Temp(t) 

thickness(t) 


Fiaure 2. Schematic of biophysical model oj sou near mansje, — — 

inmate variables drive ihe model, and the ouipu! Is Ihe p-wslcal slate o/ 

the modeled soil layers . 

decomposition vs. production to changes In temperature and mois- 
ture. Whether any ecosystem produces a net release or uptake of 
carbon to or from the atmosphere In a given year depends on th 
ratio of primary production to decomposition. In the Arctic, s s a 
quotient of two relatively small numbers, so that the potential for a 

sizeable change in this ratio is high. , . 

If inputs and effluxes of carbon varied in a similar fashion in 
response to climate-related perturbations, there would be less cause 
for concern, but most current evidence suggests that this is not the 
case In a recent review of soil respiration (Raich and Schlesinger n 
preparation), the authors point out that the vast majonty of studi^ 

show soil C0 2 efflux to be related to temper^ 

fashion with a mean Q 10 of 2.4. Nadelhoffer et al. (1991) found 

essentially no change in soil respiration with temperatures between 
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2° and 10°C, but a large and rapid increase above 10°. Temperature 
effects on production of Arctic vegetation are poorly understood, but 
most existing evidence suggests a more linear relationship. If this is 
the case, warming in the northern latitudes could actually result in 
carbon accumulation in the soils, thus creating a negative feedback 
to further warming. The accumulation of carbon in tundra during 
warmer, interglacial periods lends some credence to this hypothesis. 
This potential sink of carbon could only exist, however, over a lim- 
ited range; if temperatures increase enough to push the exponen- 
tially responding decomposition processes beyond those of produc- 
tion, the balance could shift to a net efflux, thereby generating a 
positive feedback leading to even higher temperatures (Townsend et 
al., in preparation). The actual responses of decomposition and pro- 
duction in Arctic environments to temperature changes and the cur- 
rent rates of these processes must be identified before accurate pre- 
dictions of future dynamics are possible. 

Temperature represents only one of many factors that can influ- 
ence carbon balance. A similar analysis may be applied to moisture, 
whose controls and likely changes in a doubled-C0 2 climate are 
even less well understood. In general, both decomposition and pro- 
duction often increase in an exponential fashion with greater mois- 
ture availability over certain ranges, but excess moisture will render 
soils anaerobic and retard rates of production and decomposition. 
Due to the nonlinear nature of these processes, the responses and 
current conditions need to be worked out. 

These challenges, however, do not preclude the design and pre- 
liminary application of a simulation model for carbon dynamics. Fig- 
ure 3 is a diagram of the model’s general structure, with total 
ecosystem carbon being divided into three vegetation compartments 
and three soil compartments. At present. Century calculates pro- 
duction as a function of temperature and moisture; we believe Arctic 
vegetation may also require an irradiance parameter. On shorter 
time scales and in the tundra, woody biomass should have little 
impact on the climate-induced responses; thus it is represented as 
an isolated pool. Both root and aboveground biomass create residue 
organic matter, which in tum feeds into the available soil carbon 
pool. The third soil pool, permafrost, is completely recalcitrant, but 
can feed into the other soil pools upon thawing. The details of this 
exchange will be determined by the biophysical model outlined ear- 
lier. Estimates of allocation between above- and below-ground mate- 
rial can be made from a variety of studies (cf. Chapin et al 1986a 
1986b; Chapin and Shaver, 1988), as can estimates of lignin N 
ratios. Since most of the tundra is highly organic, soil texture will 
not play as important a role as it does in other ecosystems. Signifi- 
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C0 2 



Figure 3. Conceptual model diagram for sol l carbon, C0 2 , and CH 4 cycling in 
Arctic ecosystems . 


cant transport of organic acids through peat layers and into streams 
and groundwater could occur. 

Waterlogged conditions may slow decomposition and the subse- 
quent release of C0 2 . but they do allow methanogenesis to occur; a 
warmer, wetter Arctic could also result in a positive feedback to fur- 
ther warming via accentuated CH 4 release from anaerobic peatlands 
and bogs. Work is currently under way to expand the Century struc- 
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ure to simulate CH 4 fluxes. Briefly, carbon in the soil and residue 
p oo s may become available for methanogenesis when anaerobic 
cond turns arise. The proportion of this carbon that is converted to 
CH 4 is a function of soil temperature and redox status. Redox itself 
s a function of moisture and microbial activity, so that microbial 
utihzation of oxygen in the previous time step will affect the redox 
stete of the present time step. Oxidation of CH 4 produced to C0 2 is 
a so determined by redox state, as well as by the speed at which 
CH 4 can escape to the atmosphere. A greater percentage of CH 4 pro- 
duced will be released in vegetated areas, because plants provide 
ective conduits through the water/soil layer (Conrad 1989- 
Schut, et al„ 1989a. 1989b,. Where vegetation does no, supersede 
s anding water, transport from source areas to the atmosphere is 
slower relying on diffusion, ebullition, and wind -generated turbu- 

Ih nCe iI°K C ! Cape - F ° r thIS reason ’ ve g etated ^d open-water areas 
° U d ’jf treated differently in calculations of CH 4 flux. Finally, 
H 4 is allowed to enter the soil/water environment from the atmos- 
phere as a function of the concentration gradient. There is some 
suggestion (Steudler et al„ 1989; Mosier et al„ 1991) that CH 
uptahe by soils may be closely tied to N status, so that it may even 

U , a y K e u?° SSible tG tIC ^ flux to the Mineralization rates calcu- 
lated by the central model. 


Regional-Scale Modeling 

iq«?i at S! he f terogeneIt y can be extreme in Arctic regions (Miller, 
1982). Therefore, extension of the model to regional scales will be 
challenging As there is some evidence that Arctic vegetation can be 
lghtly correlated with water status, some knowledge of topography 
(e.g. via digital elevation models) might also allow estimates of the 
wgetation parameters. This relationship is by no means certain; 
Miller (1982) states that vegetation zones in central Alaska are 
largely determined by nutrient availability rather than water status. 
The best approach to spatial extrapolation of the local model may 

^ °/ ^ Ct aL (1989) ’ Which uses Probability distributions 
and Monte Carlo sampling techniques to generate expected values 
roughout large regions to provide inputs to the ABISKO II (Bunnel 

and8COUllar ' 19751 modeL Estimates of the probability functions 
r the model inputs might be possible from a combination of 
regiona vegetation maps, AVHRR data, and thematic mapper data 
bases. Using a probability function approach to calculate expected 
values in space may be critical in the Arctic, as many of the impor- 
tant dynamics respond in nonlinear fashion, and conventional aver- 
agmg-in- space approaches could produce greatly erroneous results. 
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Critical Field Measurements 

Complete data sets needed for model parameterization and valida- 
UonTe scLce but some Information is available from sites in Mm- 
Alaska Scandinavia, and the Hudson's Bay region of 
Canada. Efforts to study carbon dynamics in northern 
increasing (e.g.. Svensson and Rosswall. 1984. apm d 

iQRfih- Sebacher et al., 1986; Cnll et al., 1988, Whalen ana 

Reeburgh, 1988; Nadelhoffer et al.. 1991; and 

tics and Space Administration's Global Troposphere 
Experiment/Arctic Boundary Layer-3 experiments). NASA, the 
NaUonaJ 1 Environmental Research Centre (United Kingdom), and the 
CanadUin Institute for Research in Atmospheric Chemistry are cur- 
rCTitb^^annhig large-scale field campaigns in the boreal arid tundra 
regimis during the next few years. Nevertheless, we will have to 
Sow much more to test and direct regional-scale simulation mod- 
els A list of needed information follows. 

. Production, decomposition, and 0O 2 and CH, ifl- 

coupled with soil moisture and temperature over at least g 
ing season, including freeze and thaw periods, and pre era y 
a full year in sites that span the general range of Arctic ecosys e 
types. This is the most critically needed information. 

. Time series of winter snowpack. snowfall, air temp«atu£. 
the effects of topography on snowpack. Some such data exls 

(e.g.. Woo et al., 1983). 

. Depth profiles of organic carbon to the permafrost flayer mrougl h- 
out the Arctic. Again, limited data are available, and Doolittle 
“ , ggofsuggeS that ground-penetrating radar may provide a 
I, d and accurate means for determining the depth to frozen soil 
and thus the thickness of the active layer. 

. rp. r hon-13 carbon- 14, and deuterium isotope values for gas 
fluxes. The deuterium isotopic ratio can be used to determine 
whether methane was produced from acetate or car on ioxi , 
the carbon isotopes allow differentiation among various types o 


substrate. 

. Redox profiles (pH, pe) for various soil moisture conditions. One 
set of such data is available from northern Alaska. 

. Soil texture (where applicable) and soil hydraulic and thermal 
properties for both organic and mineral soils. 

. Root-shoot ratios and the lignim.N ratios for each vegetation com- 
ponent in different Arctic regions. 
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